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Abstract 

In graphene, the valleys represent spin-like quantities and can act as a physical resource in 
valley-based electronics to novel quantum computation schemes. Here we demonstrate a direct 
route to tune and read the valley quantum states of disordered graphene by measuring the meso- 
scopic conductance fluctuations. We show that the conductance fluctuations in graphene at low 
temperatures are reduced by a factor of four when valley triplet states are gapped in the presence 
of short range potential scatterers at high carrier densities. We also show that this implies a gate 
tunable universal symmetry class which outlines a fundamental feature arising from graphene's 
unique crystal structure. 
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Quantum interference of electrons (or holes) causes the electrical conductance G of a 
disordered metal to fluctuate aperiodically yet reproducibly with Fermi energy, magnetic 
field or disorder configuration lM2j. When the sample size is smaller than the phase co- 
herence length (L^), the conductance fluctuates with a universal magnitude ~ e 2 /h, the 
quantum of conductance, irrespective of material properties, device geometry or dimension- 



ality. In conventional mesoscopic conductors, such as thin metal films [lj or weakly localized 
semiconductor devices [4], both universal conductance fluctuations (UCF) and weak local- 
ization effects are well-understood on the common platform of coherent backscattering of 
quasiparticle wave functions [2J. 

The scenario is considerably more complex in graphene due to the existence of two de- 
generate valleys [5j . The hexagonal lattice structure of graphene contains two basis atoms 
in its unit cell, which makes the conduction and valence bands of graphene cross the Fermi 
level at two inequivalent K and K' points on the opposite corners of the hexagonal Brillouin 
zone, leading to two valleys Consequently, quantum correction to the conductivity 

in graphene is determined not by the inelastic processes alone, but also by elastic collision 
events that involve the intervalley and intravalley scattering 8HlO]. 



TABLE I: Details of the single layer graphene devices 



Device 


Device area (LxW) 1 


Dirac point 2 


electron mobility 3 


hole mobility 3 


Dev I 


2.1 x 2.9 


-5.5 


12,000 


10,000 


Dev II 


2.5 x 3.6 


-8 


4,500 


5,700 


Dev III 


3.3 x 5.5 


40 


2,300 


2,800 



1 both dimensions in /jm 

2 in Volt 

3 in cm 2 /Vs at 4 K 



Although the signature of UCF in low temperature electrical transport in mesoscale 



llMl4|. a quantitative influence of 



graphene has been reported in numerous experiments 
valleys on UCF has never been observed. Numerical investigations suggest that the UCF 
magnitude in graphene should be stronger than in normal metal, raising doubts whether its 
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magnitude is "universal" at all [15j. Moreover, UCF is sensitive to the universal symmetry 
class of a disordered system which in case of graphene is determined by time reversal sym- 



metry (TRS) as well as valley degeneracy [8|, |l6|, Il7j . Thus UCF may reveal, for example, 
whether graphene belongs to the symplectic symmetry class in the absence of intervalley 
scattering, or whether time reversal symmetry is broken by ripples due to substrate rough- 
ness [18j, or edge magnetism 19 1 etc.. 

In this letter we present the first direct manifestation of the valley coherent states on 
UCF in monolayer disordered graphene. The key result is the suppression of UCF mag- 
nitude within individual phase coherent boxes by a factor of four as the carrier density is 
increased from close to the Dirac point, where both valley singlet and triplets contribute to 
the fluctuations, to the high electron or hole density regime, where short range potential 
fluctuations allow only valley singlets to survive. We find the factor of four suppression to 
be universal, irrespective of the geometry, carrier mobility or temperature, indicating it to 
be a unique and fundamental property of mesoscopic graphene. 

Graphene devices were prepared by standard micro-mechanical exfoliation of natural 
graphite on 300 nm Si02/Si ++ wafer surface (se e the SEM image of a typical graphene device 
in the inset of Fig. la and our earlier works 20j for more detail on device fabrication). Here 
we present detailed studies on three monolayer graphene devices with varying disorder (see 
Table I for details). The conventional magnetotransport was studied first, where we used 
the method described in Ref. to calculate the average magnetoconductance (MC) within 
small gate voltage (Vbg) windows of 4 V for each transverse magnetic field (B). In Fig. lb 
we show the quantum correction to conductivity Aa(B) = cr(B) — cr(0), as a function of B 
in Dev I at three values of Vbg which are identified on the resistance, R — Vbg trace for 
the same device in Fig. la. Measurements were performed in two terminal configuration 
as appropriate for UCF [21], although, for fitting the B— dependence of MC we have used 
the four-terminal resistivity to eliminate the contact resistance. With decreasing density 
(moving from region 3 to region 1) MC changes its sign from positive to negative similar to 
the observations in Ref. jlO |. This transition from weak localization to weak anti-localization 
indicates that the scattering associated with short-range defects is stronger at high carrier 
density and one can tune the relative strength of the elastic scattering just by chan ging the 



carrier density or gate voltage. On fitting the traces with the theoretical expression 



221, we 



find that La. saturates to a value that is in the order of the device dimension in all cases 
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below ~ 4 K, which is about 5 — 10 times longer than the elastic mean free path 23]. 

In order to evaluate the magnitude of conductance fluctuations, we chose successive gate 
voltage windows (of equal width AVbg = 4 V) within which the average conductance does 
not vary appreciably, but we have significant fluctuations for statistically meaningful analysis 
(up to ~ 800 realizations). The conductance variation within a typical window is shown in 
Fig. lc, where the random yet reproducible fluctuations in G has an amplitude of ~ e 2 /h, 
which is the hallmark of UCF. The fluctuations weaken with increasing temperature (T), as 
illustrated in Fig. Id. Below ~ 300 mK, (d~G 2 ) becomes nearly constant as itself saturates 
to ~ device dimension [23 |. 

In order to establish the origin of the conductance fluctuations, we measured the fluctua- 
tion magnitude in a small transverse magnetic field. In Fig. 2 the variation in the normalized 
noise magnitude N(B) = (5G(B) 2 ) / (G) 2 with B is shown for Dev III at three values of Vbg 
extending from the Dirac region to high carrier densities (see inset). In all cases N(B) 
decreases by ~ factor of two, as B exceeds ~ 10 mT which is also the characteristic field 
scale associated with the quantum correction to conductivity (B^ = h/eL^ « 1 — 4 mT 
for Dev III). In the diagrammatic representation of quantum transport, this can be read- 
ily explained by the suppression of the Cooperon (particle-particle channel) contribution to 
UCF for B ^> while the diffuson (particle-hole channel) contribution, which is equal in 

n n 

magnitude, remains unaffected |2J, |24J. This result has two important implications: First, 
in the framework of the random matrix theory, the factor of two reduction in UCF can 
be understood by a S-induced lifting of TRS that drives the system from orthogonal or 
symplectic symmetry class to the unitary class jsj. Similar observation here indicates that 
at B = there is no spontaneous breaking of TRS in our graphene devices. Second, we 
conclude that the observed Vec-dependent fluctuation in G is entirely due to UCF, and 
hence the measured magnitude of (5G 2 ) can be directly used to probe the valley effects and 
symmetry class of graphene. 

Since valley degeneracy is connected to the carrier density through the nature of potential 
scattering {5, 8, 22], we subsequently measured (5G 2 ) over a wide range of carrier density in 
all three devices at T = 10 mK and B — 0. Fig. 3a shows the variation of (5G 2 ) with Vbg 
for Dev I, indicating a weak increase on both sides of the Dirac point. The variation of (5G 2 ) 
with Vbg was found to be highly device specific, with opposing trends observed in different 



devices 



231 ] . The measured (SG 2 ) however arises from the entire sample, and in order to 
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e phase coherent box (L^ x £,</,), we have 



{5G 2 )/(G) 2 = (1/N) x (5G 2 )/(G^ 



estimate the UCF magnitude (SGffy within a sing' 
employed the theorem of classical superposition 
where N = LW / L 2 is the number of phase coherent boxes in a device of length L and width 
W. Using (Gq) = a, the conductivity and (G) = Lcr/W, we get, 



(SGI) = 777 x 



(6G 2 ) 



(1) 



W L 2 ■ 

In Eq. [U both (5G 2 ) and L^> are Vbg~ dependent quantities, and we evaluated them experi- 
mentally to extract the Vbg - dependence of (^G 2 ). Fig. 3b shows the dependence of L$ with 
Vbg i n Dev I over the same range of Vbg- ^4> increases with increasing density, and near the 
Dirac region it shows a minimum. This is typical for graphene and could be connected to 
dominance of electron-electron scattering with increasing densities 25]. In Fig. 3c we have 
shown the variation of (SGV) with Vbg m Dev I at B = 0, using experimentally measured 
(5G 2 ) and (from MC) at every gate voltage. While the absolute magnitude of (<5G 2 ) is 
expectedly of the order of (e 2 /h) 2 in all cases, the key aspect is the suppression of (SG 2 ) by 
a factor of ~ 4 as Vbg is varied from the Dirac point towards both high electron or hole den- 
sity regime. The characteristic scale of such a suppression seems to follow the crossover of 
linear (Coulomb scattering) to sub-linear (short-range scattering) 2f| density-dependence of 
conductivity (see Fig. 3c and 3d, and vertical guidelines). Similar reductions were observed 
in the other two devices (Dev II and Dev III) as well, shown in Fig. 4a-d. The suppression 
of noise is often asymmetric about the Dirac point, probably connected to the difference in 
mobility and disorder between the electron and hole side. 

TABLE II: Values of the symmetry parameters k, s, f3 for different universality class in case of 
graphene. 





Density 


s 


k 


P 


ks 2 

P 


Universality class 


valley-isospin 


Kramer's dej 


;eneracy 


real spin 


B = T 


n — > low 
n — > high 


2 
1 


2 
1 


2 
2 


1 
1 


4 
1 


16 

4 


GSE 
GOE 




n — > low 


2 


1 




2 


1 


2 


8 


GUE 


B = 100 mT 






















n — > high 


1 


1 




2 


1 


2 


2 


GUE 



The factor of four suppression in (5G 2 ), irrespective of mobility, geometry and other de- 
vice specific details, suggests a fundamental property of disordered graphene. In a recent 
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analytical approach, removal of valley degeneracy has been shown to cause an exact factor 
of four reduction in the mesoscopic fluctuations in graphene Q. The underlying physical 
mechanism can be schematically presented as in Fig. 4e. At lower density, the scattering is 
dominated by the long range Coulomb scattering where the effective valley-isospin rotational 
symmetry (SRS) is preserved. Consequently, the mesoscopic fluctuations receive equal con- 
tribution from each of the singlet and (three) triplet channels of diffusons and Cooperons. 
However, at higher densities the intervalley scattering dominates, and the effective SRS is 
lifted, and only the singlet diffuson and Cooperon channels contribute. In the latter case 
the graphene essentially behaves as a conventional 2D disordered metal. 

The existence of valleys is also expected to result in a nontrivial symmetry property of 
graphene Hamiltonian that has never been probed in a direct manner. At zero magnetic 
field the presence of valley isospin rotational symmetry (no intervalley scattering) makes 
graphene belong to the Gaussian symplectic ensemble (GSE), characterized by a Wigner- 
Dyson parameter f3 = 4, and doubly degenerate isospin and effective time reversed states 
(Kramer's degeneracy). Removal of valley degeneracy results in a Gaussian orthogonal 



ensemble (GOE) with /3 = 1 as well as suppression of effective Kramer's degeneracy [16 ]. 
For open mesoscopic systems, the random matrix theory connects the magnitude of UCF to 
its symmetry class simply as jsj], 

where k is the number of independent eigenvalue sequences of the transmission matrix or 
Hamiltonian, s is the eigenvalue degeneracy. As shown in the first row of Table II, factor of 
four suppression of UCF with increasing carrier density can be explained from the change 
in symmetry properties of graphene as the valley degeneracy is removed. 

Finally, in order to establish the consistency of Eq. |2l we have measured the Vbg de- 
pendence of (SGV) at a finite B (~^> -^j)-, which removes the time reversal symmetry and 
puts graphene in the Gaussian unitary ensemble (GUE, (3 = 2) at all densities. As shown 
for Dev III at different temperatures (Fig. 5), (SGV) is again suppressed by a factor of four 
at high densities irrespective of T, which can be understood either as gapping of triplet 
diffusons, or alternatively as, suppression of valley isospin degeneracy in Eq. [2] (see lower 
row of Table II). 

In conclusion, we have demonstrated for the first time that the mesoscopic conductance 



fluctuations in single layer graphene are dependent on the valley hybridized states and phys- 
ical symmetries of Hamiltonian. We have shown that the variance of universal conductance 
fluctuations can be a sensitive probe to this, which in a single phase coherent box increases 
by four times near the Dirac region as compared to high density. This sensitivity could be 
exploited in read-out schemes involving valley states in graphene. 
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FIG. 1: (color online), (a) Resistance (R) vs. gate voltage (Vbg) characteristics for a single layer 
device (Dev I) at T = 10 mK. Inset shows the schematic along with the false color SEM image 
of a typical graphene transistor, (b) Evolution of Magnetoconductance with increa sing electron 



22 



m- (c) 



density measured at T = 4.5 K. Solid lines are fit to the weak localization theory 
Conductance fluctuations as a function of gate voltage for Dev I, shown for various temperatures. 
Curves at different temperatures are shifted for clarity, (d) Variance in conductance, calculated 
from the conductance fluctuations in gate voltage range —50 V to —46 V for Dev I, plotted as a 
function of temperature. 
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Figure 2 



FIG. 2: (color online). Magnetic field dependence of the magnitude of conductance fluctuations 
for Dev III at T = 10 mK at various gate voltages: (a) 42.5 V, (b) 12.5 V, and (c) -34 V. For 
comparison, we have plotted the ratio N(B)/N(B = 0), with N(B) = (5G 2 )/(G) 2 at the magnetic 
field B. Inset shows the R — Vbg characteristics of Dev III at T = 10 mK with arrows indicating 
the three measured regions. 
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FIG. 3: (color online), (a) Variance in conductance ((5G 2 )) vs. gate voltage (Vbg) at 10 mK 
for Dev I. (b) Gate voltage dependence of and Li at T = 10 mK, extracted from the weak 
localization fits, are shown for Dev I. (c) The variance in conductance in a phase coherent box of 
area L?, (SG^), extracted from (a) and (b), is plotted with gate voltage for Dev I. The factor of 
four is indicated by the dashed lines, (d) Conductance (G) vs. gate voltage (Vbg) for the device at 
T = 10 mK. The solid line indicates the linear region in both electron and hole sides. The vertical 
dotted lines in (c) and (d) indicate the densities where the short-range scattering dominate. The 
shaded region indicates the inhomogeneous region near the Dirac point. 
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FIG. 4: (color online), (a)-(b) Conductance (G) vs. gate voltage (Vbg) at 10 mK for Dev II 
and Dev III, respectively. The solid line indicates the linear region in both electron and hole 
sides, (c)-(d) The variance from conductance fluctuations in a phase coherent box of area is 
plotted with gate voltage for Dev II and III, see details in Table I. The factor of four is indicated 
by the dashed lines. The vertical dotted lines indicate the densities beyond which the short- 
range scattering dominates, (e) Schematic describing the crossover from orthogonal to symplectic 
universality class depending on the presence of short range scattering which breaks the effective 
valley-isospin rotational symmetry with (5Gq) being the variance in conductance for a single phase 
coherent box at very high density. 
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FIG. 5: (color online), (a)-(d) Gate voltage dependence of the ratio {SGfi) / {SGq) for Dev III, 
plotted for various temperatures at B = 100 mT, where (8Gq) is the variance in the conductance 
for a single phase coherent box at Vbg = — 68 V. 
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MAGNETOCONDUCTIVITY MEASUREMENTS 

For magnetoconductivity (MC) measurements we have averaged out the fluctuations over a gate voltage window 
of AVbg = 4 V, following Ref. [H, to get the conductivity correction Act = {o(Vbg,B) — ct( Vbg, 0))av bg ■ The 
quantum correction depends not only on the dephasing rate tT , but also intervalley (t" 1 ) and intravalley scattering 
rate (r" 1 ). To fit the MC data we have used the weak localization theory for graphene, given in Ref. |2j, 

Here, F(z) — Inz + -0(0.5 + z~ v ), ip(x) is the digamma function, r^ 1 = 4eDB/h [D is the diffusion constant]. 
The diffusion constant, D — vpl/2, was determined from the mean free path I = h/2e 2 kpp, where p is the four 
terminal resistivity at that gate voltage. The corresponding scattering lengths can be calculated using the formula, 

Temperature dependence of scattering lengths 




B ( mT ) T (K) 

Fig. SI. (a) Magnetoconductivity as a function of temperature taken at an electron density n — 1 x 10 12 /cm 2 . 
Solid lines are fit to the weak localization theory (Eq.Q}. (b) The scattering length scales as a function of temperature, 
extracted from the fits. 

Fig. Sla shows the MC data at a particular density (f x 10 12 /cm 2 ) for various temperatures for Dev I and the 
corresponding fits are shown by solid black lines. With increasing temperature the inelastic phase coherence length 
(Lcf,) decreases (Fig. Sib) and saturates at lower temperatures (< 4 K) at this particular carrier density. However, 
the intervalley and intravalley scattering lengths remain almost constant with temperature as these are related to the 
elastic scattering 0, Q ■ 

Gate voltage dependence of and Li 

Fig. S2a and S2b show the gate voltage/density dependence of MC data at the base temperature of ~ 10 mK for 
Dev II and Dev III, respectively. It is clear from both the graphs that with increasing density L<p increases, whereas, 
the intervalley scattering length (Li) remains unchanged. The density dependence of the scattering lengths have been 
studied in detail in Ref. To calculate the fluctuations from a single phase coherent box we have divided the area 
between the voltage probes in small boxes of area L 2 , : schematically shown in Fig. S2c. We have used the L^ value 
for a particular density, extracted from the fits. 
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Fig. S2. Gate voltage dependence of and Li at T = 10 mK, extracted from the weak localization fits, are shown 
for (a) Dev II and (b) Dev III. (c) Schematic describing that the area (L x W) between the voltage probes can be 
divided into phase coherent boxes of area L\. 

GATE VOLTAGE DEPENDENCE OF (SG 2 ) 

Fig. S3a and S3b shows the variation of (SG 2 ) with gate voltage for Dev II and Dev III, respectively. (6G 2 ) was 
calculated using the formula (SG 2 ) = (SR 2 ) / (R) 4 , where (R) is the average four terminal resistance in that gate 
voltage window. We observe that (5G 2 ) is very much device specific and hence, is not a significant parameter to 
understand the absolute magnitude of UCF. 
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Fig. S3. Variation of (SG 2 ) with gate voltage at T = 10 mK, for (a) Dev II and (b) Dev III. 
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